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Abstract 

We analyse the spatial diversity of a multipath fading process for a finite region or curve in 
04 ' the plane. By means of the Karhunen-Loeve (KL) expansion, this diversity can be character- 

ful ised by the eigenvalue spectrum of the spatial autocorrelation kernel. This justifies to use the 

. term diversity spectrum for it. We show how the diversity spectrum can be calculated for any 

' such geometrical object and any fading statistics represented by the power azimuth spectrum 

(PAS). We give rigorous estimates for the accuracy of the numerically calculated eigenvalues. 
The numerically calculated diversity spectra provide useful hints for the optimisation of the 
geometry of an antenna array. Furthermore, for a channel coded system, they allow to evaluate 
the time interleaving depth that is necessary to exploit the diversity gain of the code. 
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q ; 1 Introduction 



Multipath fading is a major source of degradations in wireless digital communication systems 
[Q [H EH- Information bits that fall into deep fades may get lost if no measures are taken against 
this. All methods to prevent this loss of information employ some kind of diversity, i.e. they 
spread the information over a wider range of the channel so that enough reliable parts of it can 
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■^j- ■ be exploited by the transmission system. In its simplest form, diversity is employed by multiple 

reception of the same information bits at different frequencies, time slots, or antenna locations. 
More efficient ways to employ diversity is channel coding in conjunction with an appropriate time 
and frequency interleaving [T1E3I3JEI] & s applied in several systems. Multiple- input- multiple-output 
(MIMO) systems with several transmit and receive antennas are an efficient method to increase the 
channel capacity [3J[5J[7]- To achieve the maximal capacity, these antennas must be uncorrelated, 
which requires a sufficient spatial distance between them. 

x ■ . . , 

' In practice however, the information can only be spread over a finite region of the channel. There- 



fore, there are channel correlations that diminish the possible gains of all of these methods. In 
this paper we provide tools to quantify these degradations. We concentrate on space and time 
diversity aspects of the multipath channel and leave aside frequency diversity as a topic of its own. 
Space and time diversity, however, are closely related for the scenario of a receiver that is moving 
through a spatial fading pattern: If the receive antenna moves along a straight line, the time fading 
amplitude as a function of time is nothing but the fading amplitude as a function of the space 
variable, scaled by the velocity of the vehicle. Thus, time correlations are nothing but scaled space 
correlations. The diversity of such a given space or time interval can be quantified by evaluating 
the equivalent uncorrelated diversity branches of the fading process. For a time interval [ii,i2], 
these uncorrelated branches can be obtained by the Karhunen-Loeve (KL) expansion [51 03 [TO] of 
the fading process a (t) given by 

oo 

a (t) = J2 c ^n (t) , t 6 [h,t 2 ] . (1) 
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In that expansion, the basis functions ip n (t) are orthogonal and of normalised energy, and the 
coefficients c„ obey the condition 

E{ (2) 

This means that the process can be decomposed into a series of orthogonal and uncorrelated com- 
ponents, each with energy X n . We may thus call them equivalent uncorrelated diversity branches. 
Their respective strenghts are characterised by the values A„ , which are the eigenvalues of an in- 
tegral equation with kernel given by the autocorrelation function ( ACF) of a(t) [FJ [HI OH] • The 
eigenvalue spectrum {An}™, characterises the amount of diversity gain that can be achieved in 
the time interval under consideration. We will thus call {X n }^Li the diversity spectrum of the 
stochastic process corresponding to that interval. Coding with interleaving is closely related to 
diversity. The error event probabilities corresponding to the Hamming distance d of a convolu- 
tional code with perfect interleaving in a Rayleigh fading channel are asymptotically proportional 
to SNR~ d , i.e. they behave like the error probabilities of d uncorrelated diversity branches [TJ[51[3]. 
For interleaving over a finite time interval, this error probability will decrease much slower if there 
are less than d eigenvalues of significant size 

In that example where we consider time variance that is caused by vehicle motion, the diversity 
spectrum can equivalently and uniquely assigned to a corresponding spatial interval (the segment 
of a straigth line). On the other hand, a spatial interval by itself can be interpreted as the 
geometrical base of a linear antenna array. In that case the term (antenna) aperture is usually 
applied. With the continuous-limit idealisation that the straigth line is densely covered with 
sensors, the eigenvalue spectrum may then be interpreted as the multipath richness of the fading 
process (or the electromagnetic field) in that region. We may thus say that the eigenvalue analysis 
of a time interval and densely covered uniform linear array (ULA) are equivalent. 

Our goal is to calculate the diversity spectrum not only for a given spatial (or time) interval as 
discussed above, but for quite general one- or two-dimensional apertures in the plane. The two- 
dimensional aperture of a rectangle, e.g., arises naturally if a linear antenna array is mounted on 
a vehicle, perpendicular to the direction of motion. One edge of the rectangle corresponds to the 
length of the array and the perpendicular edge corresponds to the route covered by the vehicle 
during a given time interleaving depth. 

Rather than characterising the diversity of a given aperture by the complete diversity spectrum 
{A n }^L-,, it is desirable to have a single number that reflects the diversity degree and can be 
interpreted in some way as the effective number of equivalent independent diversity branches. It 
turns out that the diversity measure defined by 

is a useful quantity that allows such an interpretation. This definition is the dense-array limit of 
the L-antenna diversity measure 

that has been discussed in-depth by [H] and Q15]. For MIMO systems, one can calculate the 
diversity measures lutx and Drx at the transmitter (Tx) and receiver (Rx) site, respectively. Under 
the assumption that there are no correlations between both sites, the spectral efficiency slope Eo 
of the low power regime |14l I13| is given by 

2 11 

tt = + (5) 

^0 ^Rx ^>Tx 

Our method allows to calculate So for dense arrays. 

The major part of our work is devoted to the solution of the eigenvalue problem in the Hilbert 
space of square integrable functions over a given aperture. In |15j . this problem has been solved for 
the special case where the aperture is a disk of finite radius. The authors approximate the fading 
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process by a truncated orthogonal series expansion based on Bessel functions, and they prove trun- 
cation error bounds for the fading process. Based on that rigorous bounds, they heuristically argue 
that a good approximation of the process itself leads to a good approximation of the eigenvalues 
An . 

We adopt the method of [TS] and put it into a more general Hilbert space framework. This allows 
us to calculate diversity spectra not only for a disk, but for quite general curves or areas in the 
two-dimensional plane. Furthermore, by utilising the theory of compact operators, we are able 
to prove rigorous bounds of the approximation errors for the eigenvalues themselves. We give 
several examples of numerically calculated eigenvalue spectra and diversity measure for scenarios 
of several apertures and scattering environments with azimuth spectra (PAS) of different degrees 
of anisotropy. 

The rest of the paper is organised as follows. In Section 2, we introduce the channel model and 
discuss its statistical properties. In Section 3, we introduce the Hilbert space of square integrable 
functions on the aperture and define the the autocorrelation kernel as an operator acting on 
that Hilbert space. Numerical examples of diversity spectra and diversity measures for different 
scenarios are presented in Section 4. In Section 5, we summarise our results and draw some 
conclusions. The necessary Hilbert space framework of our work as well as the proofs of the 
theorems can be found in the Appendix. 



2 The Statistical Model for Spatial Multipath Fading 
2.1 Solutions of the Helmholtz Equation 

We consider scalar spatial multipath interference (fading) patterns that are given as scalar solutions 
4> of the two-dimensional (2D) Helmholtz equation 

A0 + fc 2 = O, (6) 

where A= d 2 / dx 2 + d 2 / dy 2 is the 2D Laplacian operator and k is the wave number. The Helmholtz 
equation (HE) is equivalent to the homogeneous wave equation for a fixed wavelength. Considering 
only the homogeneous equation means that all scatterers are outside the region of interest. Since 
we keep the wavelength fixed, we may set it equal to one which means that we measure the space 
variable x £ R 2 in units of wavelengths. With that normalisation, the wave number is simply given 
by k = 2t:. 

The prototype solution of the HE is given by 

0(x) = e>" k ~, (7) 

with |k| = k = 2n. The wave vector — k points into the direction of propagation, thus k points 
into the direction of the incoming wave. We write k and x in polar coordinates as 

k = 27 rf C ° Sa ) and x = r ( ^ ) , (8) 

Y sin a J \ smp J w 

note that 

k • x = 2irr cos (a — 0) (9) 

and obtain 

0(x) = c 3^rco S { a -f3) _ ( 1Q ) 

The angle a is the angle of arrival (AoA). 

A general solution a (x) of the HE can be obtained as a superposition of solutions of type (flU)) 
with different angles of arrival a: 

o(x) = ^ a(a)e^ rcos(Q -^da (11) 

J — 7T 
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In that equation, a (a) is a weighting factor, and a (a) da has to be interpreted as the part of the 
fading amplitude a (x) corresponding to AoAs corresponding to the directions inside the interval 
[a, a + da]. We expand the weighting factor a (a) into a Fourier series 



with coefficents 



a [a) 

n— — oo 



~ jan &(a)^. (13) 



For convenience, we define 



a n = 2na n = e ^ an a (a) da . 

J — 7T 

Inserting the Fourier series (|T12j) into Equation (|TT|) we obtain: 

o(x) = / a(a)e* 2,rrcoB ( a -^da 

J— 7T 

OO 

= ^ a n / e J ' an e 7 ' 2,rrcos(Q!- ^ ) da 

71 — — OO 

oo 

= ^ a„ / e j(Q+/3)? V 27rrcos(Q) da 

n=-oo u '^ 7r 

OO 

^ a„c j/3 " / e jan e 32 " cos(a, (ia 



71— — OO 



(14) 



J„(x) = -J_ / e JQ "e J:ccos(a) da, (15) 



The complex version of Bessel's first integral ([12], p. 690) 

1 

27Tj" 

then leads to the series expansion 

OO 

o(x)= a n ^ n j n J n (2wr) . (16) 



2.2 Properties of the Fading Random Process 

The fading amplitude a (x) is modeled as a zero-mean stochastic process. This means that a (a) 
in Equation (fTTj) is a family of random variables depending on the angle a. We impose the 
uncorrelated scattering assumption (US) [171118) which means that waves impinging from different 
directions can be regarded as uncorrelated: 

E {a (a) a* {a')} = S (a) 5 (a — a') (17) 

The quantity S (a) is called the power azimuth spectrum (PAS), and S (a) da has to be interpreted 
as the amount of average power that is impinging from directions inside the interval [a, a + da]. 
The PAS must be positive and integrable. We normalize the PAS to the total power one: 

S((p)d<p = l (18) 

From Equation (fTT|) , together with the US assumption (TIT)) , we obtain the following expression for 
the autocorrelation of the fading process: 



E{a(x)a* (x')} = I p) e -i2«r' cos( Q -p') s ^ da 

J — 7T 



(19) 
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Using Equation we note that 

r cos (a — P) — r' cos (a — /?') = k • (x — x') . 
Thus, E {a (x) a* (x')} depends only on the difference variable, x — x , and we may write 

E{a(x)a* (x')} = p (x - x') (20) 
with the spatial autocorrelation function p (x) given by 

p(x)= / c J ' 27rrcos(Q -' 3) t S(a)dQ ! . (21) 



From Equation l|18p we obtain p (0) = S (ip) dtp = 1 . Thus, the autocorrelation function has 
the property 

|p(x)|<l. (22) 

The PAS has a Fourier expansion 

oo 

S(a)= S « eJan ( 23 ) 

71 — ~ OO 

with coefficients given by 

e-*»S(a)^. (24) 

Z7T 

We assume that S (a) is a piecewise continuous, bounded function over the interval [— tt, it]. Then 
the Fourier series (|23p for S (a) converges almost everywhere. For convenience, we define 

s n = 2ns n = f e~i an S {a) da . (25) 

In the following, these coefficients are assumed to be known for the PAS under consideration. 
Because S (a) must be real, 

s„ = s*_ n (26) 

holds. Furthermore we have 

|« n |<*0 = l- (27) 

A series expansion for p(x) can be obtained in a similar way as before for a(x): We insert the 
Fourier expansion (|23| into Equation (|2T|l and use the integral representation (|T5|) . This yields 



p(x)= ^ S„c^ n j"J„ (2vrr) . (28) 



n— — oo 



The US property (fT7)) has consequences also for the correlations between the (random) Fourier 
coefficients a n of a (a). The autocorrelation matrix R of the discrete random process of the 
coefficients a n — 2na n is defined by its elements 

R mn = E {a m a* n } . (29) 

We insert (|25p into this equation, use the US assumption (|17|) . compare with (|23|) to obtain the 
following expression: 

Rmn — Srn—n (^^) 

The (infinite) autocorrelation matrix R can be written compactly by defining a as a column vector 
of infinite length build from the coefficients d n . Then the autocorrelation matrix can be written as 

R = E{aa t } . (31) 

The dagger f denotes the Hermitian conjucate, i.e. the complex conjugate transpose of the vector. 
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2.3 Time Correlations and Doppler Spectrum 

A single-antenna receiver mounted on a vehicle experiences time variant fading given by a(x(t)), 
where x (t) is the antenna position at time t. The time autocorrelation of the fading is given by 

E {a (x (*)) a* (x (*'))} = p (x (t) - x (f )) . (32) 

For for x (f) = vt, i.e. a vehicle moving with constant velocity vector v, 

E {a (x (t)) a* (x (i'))} = Ph™ (t - (33) 

holds with a time autocorrelation function defined by 

Ptirne (*) = P (vi) . (34) 

Consider a motion with constant velocity v max = |v| along the x-axis. According to Equation 
([8]), we have to replace rcos(a — f3) — v max t cos (a) in Equation (|2ip which leads to the integral 
expression 

Pume (t) = c^— tcos ^S (a) da . (35) 



We split up this integral according to J" da = J da + JT da, substitute a by the Doppler 
frequencjQ 

v = v max cos (a) , (36) 

and obtain 

/"max 
e^^S Dopp i er {v) dv (37) 

with the Doppler spectrum 

5 (a (!/))+ $(-(*(,/)) 

ODoppler \V) — . = . (60) 

V V max ~ V 

Here we have defined 

a (is) = arccos ( — - — ) . (39) 

It is worth to note that Equation (|38[) reflects the mirror symmetry of the Doppler shift: Two 
signals impinging from angles a and —a lead to the same value of v. Measurements along a 
straight line cannot distinguish between signals from left and right. 

Setting x = (y max t, 0) T in Equation ([28)1 yields the series expansion 

oo 

Ptirne (t) = ^ S„j"J„ {2-KV max t) . (40) 

71— — OO 

The special case of isotropic scattering with PAS 

S(a) = -L (41) 

is a popular fading channel model. The corresponding so-called Jakes [19| Doppler spectrum can 
be obtained from Equation ([3"5[) as 

SDoppler {v) = , = • (42) 

V max 

The corresponding time autocorrelation function is given by 

Ptirne (t) = Jo (I^Vmaxt) (43) 

because in the case of the isotropic PAS the series (|4TJ[) reduces to the term with n — 0. We 
emphasize that our investigations are not restricted to the isotropic scattering case. 

^^Note that, because the wavelength is set equal to one, the maximal Doppler frequency and velocity have the 
same value u max = |v|. 
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3 The Autocorrelation Operator and its Diversity Spectrum 



The goal of this section is to calculate the diversity spectrum of the KL expansion, i.e. the 
eigenvalues of the integral kernel given by the autocorrelation function of the stochastic process 
under consideration. 

Consider for example a stochastic process a (t) with time variable t restricted to a time interval 
[t±, t?\. The kernel R (t, t') = E {a (t) a* (t 1 )} defines a linear operator on the Hilbert space of square 
integrable function over the interval [ii,^]- The task is to calculate its diversity spectrum, which 
depends on the interval and on the statistics of a (t) . The KL theory for such a time interval can 
be found in Sec. 6-4 and Appendix 2 of [9]. In the following, we shall need the KL theory for a 
more general aperture given by a ID curve or a 2D region in the plane. We shall define the Hilbert 
space of square integrable functions over such a general aperture and show a method to calculate 
the corresponding diversity spectrum. 

3.1 Construction of the Hilbert Space 

Consider the following geometrical configurations called apertures and denoted by A: 

1. Continuous and piecewise smooth curves of finite length in the plane. Examples are the circle 
i={(i,i/)el 2 : x 2 + y 2 = 1 } and the interval A = {(x,y) el 2 : < x < 1 and y = 0} . 

2. Arbitrary finite and closed (topologically speaking: compact) two-dimensional regions in 
the plane. Examples are the disk A — {(x, y) G R 2 : x 2 + y 2 < l} and the square A = 
{(x,y) el 2 : x,y<l/2}. 

Both cases can be be related to an idealised antenna array with infinitely dense antenna packing on 
the curve or region, respectively. We shall analyse the multipath richness of the stochastic process 
a (x) resticted to that aperture. This has to be done by solving the eigenvalue problem of the 
autocorrelation integral kernel which is related to the KL expansion on that aperture. The first 
task is to construct the Hilbert space T~L = C 2 (A, d/x) of square-integrable functions <fi, ip, ... that 
live on that aperture A: 

[ \4>\ 2 &LL < oo (44) 

J A 

Here /j, is the (normalised) Lebesgue measure on A. The scalar product of two vectors <f>, ip G T-L is 
given by 

(4>,1>) = I (45) 

J A 

and the norm of cj> G H is defined by 

W\\ = y/fafi. (46) 

To obtain the Karhunen-Loeve (KL) expansion, we need to define the autocorrelation operator R 
on that Hilbert space. For convenience, we have normalised the measure to 

H(A)± f d/x = l. (47) 
J A 

As we shall see, this has the consequence that the trace tr (R) (i.e. the sum of eigenvalues) of the 
operator R is normalised to one. 

Some facts about Hilbert space operators and definitions of their different norms are summarised 
in Appendix A. 
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3.1.1 The Hilbert Space for a One-Dimensional Curve 



Consider a finite-length, continuous, and piecewise smooth curve A in the plane. For a given 
parametrisation r i— > x (t) the curve is defined by 

A = |x(r) S E 2 : ri < r < r 2 } , 

and .A is just the image of an interval [ti , r 2 ] . Thus, the Hilbert space of square integrable functions 
over A is isomorphic to the Hilbert space of square integrable functions over that interval. In the 
following, this Hilbert space will be constructed explicitely. The line integral of a scalar function 
T over a curve A is defined as 

/ Td£ = f "t(x(t))|x(t)|<It. (48) 

J A Jn 

In that equation, x (r) is a piecewise differentiable paramatrisation of the curve A, and x (r) stands 
for the derivative with respect to the parameter r S [ti,T2]. The integral (|48p is independent of 
the choice of the parametrization. The length of the curve can be expressed by the curve integral 
over the scalar function T = 1 as 

\A\ = I d£. (49) 

J A 

A special and frequently preferred parameter is the (running) length A of the curve. In that case 
|x(A)| = 1 holds, and Equation (|4"B")) becomes 

r\A\ 

Td£ = / T(x(A))dA (50) 



I A JO 

For this paper, we will chose r to be the normalised length 

T "W (51) 

For this so-called priviledged parameter r we obtain (by using a simple substitution) 

In that equation, we have written T (r) = T (x (t)) = T (x (A (r))) to simplify the notation. We 
shall use this loose notation whenever it does not give rise to any confusion. Equation (|52[) means 
that we can express every average over the curve by an average over the unit interval. 

Bei chosing the priveledged parameter, the Hilbert space becomes isomorhic to the Hilbert space 
over the unit interval: 

U = C?{Aii) ^£ 2 ([0,l],dr) (53) 
Here = denotes isomorphism, the measure fi is given by 

d£ 

dM = R , (54) 
and the scalar product (|45[) can be written in explicit form as 

(0» = / 0*(r)^(r)dT. (55) 
Jo 
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Example 1: The Hilbert space of functions on a circle We consider a circle of radius r 
defined by 

i={xel 2 : |x| = r} . (56) 
The priviledged parametrisation is given by 

x (t) = r (cos (2itt) , sin (2ttt)) T , (57) 

and |.A| = 2ttt. The scalar product 

/ rh*ihr\P = I 

2irr 



V) = -^ / <f>*^d£ = / 0*(r)V(r)dr. 



.4 Jo 

is just the standard scalar product on the unit circle. We note that the natural orthonormaJl basis 
on the unit circle is the Fourier basis {(in}^,^ with 

u n (r) = e^ Tn . (59) 

Assuming that each realisation of the fading is an element of the Hilbert space, we can expand it 
with respect to this basis 

oo 

a (r) = ^2 a n u n (t) (60) 

n— — oo 

with coefficients given by 

a„ = (u n ,a) . (61) 
This is the familiar formula for the Fourier coefficients: 



<(r)o(r)dT (62) 

= f C - f27TTn a(T)dT (63) 
Jo 

3.1.2 The Hilbert Space for a Two-Dimensional Region 

We now consider consider a closed two-dimensional region A C M 2 of finite area \ A\. The functions 
4> (x) that satisfy the condition 

/ |0(x)| 2 d 2 x =< oo 

J A 

build the Hilbert space 

U = C 2 (A, At) , (64) 



with the measure \i given by 



d At = ^. (65) 



The scalar product of two vector <f> , tp 6 TL is defined by 

(<t>,1>) = I 0*(x)V(x)d M (x) (66) 
J A 

<t>* (x) V (x) d 2 x . (67) 



1 

JA\ 



A 



2 The basis {u n } is called orthonormal if {u m ,u n ) = 5 mn holds. 
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Example 2: The Hilbert space of functions on a disk |15| We consider a close disk of 
radius r\, i.e. 

A = {x e R 2 : |x| < n} . (68) 
It has the area \A\ — irr 2 . The scalar product is defined by 

= 7^7 / 0(x)*^(x)d 2 x. (69) 



\A Jc 

Introducing polar coordinates r and j3 according to Equation ([5]) yields 

d 2 x d/3 rdr 
~\A\ = ^rf' 

A we shall see later, an orthonormal basis for this Hilbert space on a disk is given by 

Un ( x ) = Jn(27T|x|) _ ptf(x)w _ 

3.2 The Autocorrelation Operator 

We now define the autocorrelation operator R on the one- or two-dimensional aperture A. 

One-dimensional aperture: We define the operator R acting on % = C 2 ([0, 1]) by its integral 
kernel 

i?(r,r') = p(x(r)-x(r')) (70) 

with p (x) according to Equations ([20]) and (|2ip and x (r) being the priviledged parametrisation 
of the curve. Using loose notation a (r) = a (x (t)) in Equation (|20p we may write 

R(r ) T')=E{a(r)a*(r')} . (71) 

The action of the autocorrelation operator R on a vector ip £ H is given by 

(RV)(r)= / R{t,t')^{t')At' . (72) 



We note that the kernel is not a convolution kernel (i.e. a kernel that depends only on the difference 
variable r — r') except for the special case where A is a straight line. 

Two-dimensional aperture: We define the operator R acting on T~L — C 2 (A, fi) by its integral 
kernel 

i?(x,x') =p(x-x') , (73) 

defined for x, x' G A. Again, p (x) is given by Equations (|2T)]) and ([2"TT) . The action of the 
autocorrelation operator R on a vector ip £ % is given by 



(RVO (r) = l^(x,x')^(x')d 2 x. 



(74) 



In the following treatment, we will no further distinguish between both cases. One can show by 
direct calculation from Equation ([72"]) or (|74p that 

(0,R^) =E{(0,a)(a,^)} . (75) 

holds. The values of these so-called matrix elements (0, R^) uniquely define the linear operator 
R [2D]. Thus, Equation ([75)1 can be interpreted as an alternative definition of the operator that is 
valid for both classes of apertures. 

We list some properties of the autocorrelation operator R: 
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Continuity, integrability, and compactness: From Lebesgue's dominated convergence the- 
orem and the absolute integrability of S (a) it follows that p (x) given by integral expression (f2~Tj) 
is a continuous function on A. Thus, the kernel R (r, r') (or R (x, x'), respectively) is a continuous 
function in both variables. Because any continuous function on a closed, finite set is integrable 
(and square-integrable, too) the kernel is in £ 2 (A, dp) x C 2 (A, dp,). This property ensures that 
R is a Hilbert-Schmidt operator, which also guaranties that it is compact operator (see Theorems 
VI 22 and VI 23 in [2D] ) • Because \p (x)| < 1 it can further be shown that the operator norm of R 
is bounded by ||R|| < 1. 



Symmetry and positive semi-definiteness: Because the operator R is bounded, it is defined 
on the whole Hilbert space. Because of Equation (fTTj) . the kernel is symmetric 

R(t,t') = R*(t' iT ) , (76) 

and the corresponding operator R is self-adjoint|f|: 

R = R f (77) 

The same hold with i?(r, r') replaced by i? (x, x'). From Equation (|75|) . we conclude that R is 
positive semi-definite, i.e. 

(ip, Rip) > 

holds for every ip g H. We also say that the corresponding kernel i?(-, ■) is positiv semidefinite. 



Eigenvalues and eigenvectors of R: From the properties of R stated above, the following 
properties of its eigenvalues can be infered ([20], Chap. VI): Since R is self-adjoint, there is a 
complete set of eigenvectors of R, {^i}^, with 

(<Pi,<Pk) = Sik (78) 

that build a basis of Ji. Because R is compact, the corresponding eigenvalues eig (R) = {Aj}^ 
are of finite multiplicity and have the property 

lim Aj = . 

i— >oo 

Because the operator is positive semidefinite, we know that A^ > for all i. 



Mercer's Theorem: Because the kernel R(t,t') is continuous, symmetric, and positive semi- 
definite, according to Mercer's Theorem ([H], §5.4; [52], § 98) the following series expansion for 
the kernel holds: 

oo 

R{T,T') = Y,*m(T)y*{T') (79) 

i=l 

The theorem states that the infinite series converges uniformly in r and r 1 . For the kernel R (x, x'), 
the analogous expansion hold with r, r' replaced by x, x'. 



The trace of the operator: From Mercer's theorem we conclude 

/ R(T,T)dr = ]TA, / ^(r)^(r)dr 
Jo i=1 

„1 oo 

/ i?(r,r)dr = V A, (80) 
Jo i= i 

Due to our normalisation p (0) = 1 we have R (r, r) = p (0) = 1 and, thus, 

oo 

tr(R) = 5^Ai = l. (81) 

i=l 

The same holds for r replaced by x. 



3 See Appendix A for the definition of self-adjointness. 
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Karhunen-Loeve Theorem: Let {<£n}^Li be the basis of eigenvectors of R and eig (R) = 

{-^nl^Li be the corresponding eigenvalues defined by 

R^n = A„v?„ . (82) 

Define the (random) coefficients 

Cn=((p n ,a)= ifl (t) a (r) dr . (83) 
Jo 

One can readily concluse from Equation (|82D that these coefficients must be uncorrelated: 

E {c m c* } = S mn X n . (84) 

Conversely, one can see that the basis of eigenvectors (p n of R is the only one such that the 
coefficients (|53")) fulfil the condition The KL theorem states that the expansion 

oo 

a (t) = c„<^„ (r) (85) 

n=l 

converges as a limit in mean (l.i.m.) and that the convergence is uniform in r. The same holds for 
the variable r replaced by x. ft is worth to note that the proof of convergence (see [HIS]) is based 
on Mercer's theorem which requires the continuity of the kernel. 



Diversity Measure: We refer to Appendix A for the definition of the norms II A| 



and 



||A|| 2 of an operator A. By using these norms, the diversity measure given by Equation ([3]) can 
be re-defined as 



w(R) 



i. e. as the ratio between the trace norm 



R 



R 



i2 

,2 • 

I2 



(86) 



R 



(87) 



and the Hilbert-Schmidt norm 



R 



\ 



of the operator R. Because of hierarchy of norm given by Inequality (|152p . ui (R) > 1 always holds. 

Equation (|86p is a suitable definition of the diversity measure for an integral kernel as studied in 
this paper. The equivalent for a finite-dimensional L x L correlation matrix R it has been defined 
in [12, 13 j. In the next section we shall show numerically that for L antennas distributed uniformly 
over a given aperture, the calculated diversity measure approaches that of a dense array in the 
limit L — > 00. 

The value of the quantity oj can be understood by considering two extreme cases for a finite number 
of antennas |12[ 1 1 3] : Firstly, for L uncorrelated antennas, R is equal to the identity matrix I with 
||l||^ = L 2 and ||l||_ = L, which yields u = L. Secondly, for L antennas with 100% correlation, 

Rik = 1 holds for all elements Rik of R. This means that ||R||j = L 2 and 1 1 1 1 ^ = i 2 , which 
yields u> = 1. Thus, the diversity measure uj (R) can be interpreted as the number of equivalent 
uncorrelated antennas. Note, however, that uj (R) is typically not an integer number. 

For the continuous case, due to the trace normalisation (|81|) . we have simply 



w(R) 



> 1 



R 



(89) 
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3.3 Discrete Representation of R 



To calculate the diversity spectrum eig (R) = {A n }„ =1 , one must solve the eigenvalue problem 
given by Equation (|82[) . Expressing R by its integral kernel yields integral equation 

1 

R(T,T')^ n { T ')dT' =X n ip n {r) (90) 



n 



(or the corresponding expression in the variable x). This may be solved by numerical quadrature 
methods. However, this approach is numerically cumbersome, and the accuracy is difficult to 
control. An alternative approach introduced by |15j is to switch from the continuous representation 
given by the integral kernel R(-, ■) to a discrete representation corresponding to an orthonormal 
basis {u n } of %. Then the equivalent matrisQ R with matrix elements given by 

Rmn = (u m ,B.u n ) ■ (91) 

in this new representation has the same eigenvalues. These initially unknown matrix elements have 
to be obtained from the known matrix elements R mn defined by Equation (|30[) . From Equation 
([7SD we find that 

Rmn = E {a m a* n } (92) 



holds, where 

is the coefficient of the expansion 



(«n, a) (93) 



E 



a n u„ (94) 



of the fading process with respect to this basis. We define the column vector a of infinite length 
that is build from the coefficients a n . The autocorrelation matrix that has R mn as its elements is 
then given by 

R = Ejaat} . (95) 

The eigenvalue problem can this way been reduced to a discrete (but infinite-dimensional) one. 
Before soving it, the first task is to define the appropriate basis {u n } and then relate the matrix 
R to the matrix R with elements given by Equation 1)3011 . 

To do this, we expand the fading process a according to Equation (|16p by 

oo 

a = ^2 a n v n , (96) 

71= — OO 

where v n (x) is defined by 

u n (x) =e^« n j n J„(27r|x|) . (97) 

The domain of the functions a (x) and v„ (x) is the array A, which is either a two-dimensional 
region or a one-dimensional curve. For a curve with priviledged parameter r we may write 

v n (r)=e i ^ n j n J n (2Trr(T)) (98) 

where /3 (r) = /3(x(r)) and r (r) = |x(t)| are the polar coordinates of the parametrisation. In 
general, the set {v n } is not orthogonal. For the special case where A is a interval on the x- 
axis, the functions v n (x) — j n 3 n (2ttx) do not even form a basis. This is because of the Bessel 
function property J„ (2irx) = (— 1)" J_„ (2ttx), they are not linearly independent. But in any case, 
Equation ()96p holds for every realisation of the fading process restricted on the array. Each v n can 
be expanded into the orthonormal basis {u n } by 

fro ^ ^ T^mri^m (99) 



4 The slightly lose notation to use the same symbol for the operator and its equivalent matrix representation does 
not cause problems, and it is quite usual e.g. in quantum mechanics. 
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with coefficient^] 

T mn = (u m ,v n ) . (100) 

Expanding 

CL rn (w m , ft) 

yields 

a m = 2J T mn a n . (101) 

Let T be the matrix with elements T mn . Then the above equation can be written compactly as 

a = Ta. (102) 

We insert this into the autocorrelation matrix 

R = E{aa f } 

= E|Ta(Ta) f | 
= TEjaa^T^ 

and obtain 

R = TRT t . (103) 
Making use of the property eig (AB) = cig (BA) (see [25], p. 53), we obtain 

cig(R) = cig (TRTt) 
= eig (VTR 

i.e. 

eig(R) =eig(GR) (104) 
with the Gram matrix G that is defined by 

G = T+T. (105) 

It has the matrix elements 

G m n = (Vm,V n ) . (106) 

Since G and R are known matrices, the eigenvalue problem for R can be solved by utilising 
Equation (|104j) . Before we show how this can be done with a controlled approximation, we discuss 
the two examples from above for which the transform turns out to be quite simple. 

Example 1 (ctd.): Consider the circle (f5B)) of radius r with parametrisation given by (|57[) . 
According to Equation p^|) , the fading process a (r) can be written as 

oo 

o(t)= Yl a n j n J n (27tr)^ Tn . (107) 



n— — oo 



5 In case that {d„} is a basis, Equation HIPP! defines an basis transform operator T by 

T : u n i-* v„ = Tun ■ 

It has matrix elements 
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This is an expansion 

oo 

a ( T )= 

(r) (108) 

n— — oo 

with respect to the basis 

v n (T)=j n J n (2Trr)e i2 * Tn . (109) 

This basis is already orthogonal. Is can be normalised by dividing it by a normalisation factor £„ 
that fullfills the requirement 

l£«i = lkll ■ 

We define 

f„ = i n J„ (27rr) . (110) 
We have included the phase factor j n because the basis {u n } n a = _ 00 with 

Vn (t) = S,„U n (t) (111) 

is now just the standard orthonormal Fourier basis u n (r) = c? 2 ' nTn . 

Inserting into the expansion (|108[) and comparing with Equation (|60p , we find that the correspond- 
ing coefficients are related by 

a n = i n a n ■ (112) 
Thus, for this simple example, the transformation matrix T is diagonal. The matrix elements 
Rmn = E {a m <} and R mn = E {a m a* n } are related by 

Example 2 (ctd.) [15] : Consider the circle (|6"8")l of radius r% . The fading process a (x) can 
be written as 

oo 

0(X)= anj"Jn(27T|x|)^W", (114) 

TL— — OQ 

where j3 (x) is the polar angle of the vector x. This is an expansion 

oo 

a(x)= Y a n v n (x) (115) 

71— — OO 

with respect to the basis 

u n (x) =rj„(27r|x|)c^( x )". (116) 

This basis is already orthogonal. Is can be normalised by dividing it by a normalisation factor £ n 
that fullfills the requirement 

\U 2 = \\vnf= I ^J^27r|x|)=jT^J*(2*r). 



.4 



We define 



^ =i iL iAr Ji{27Tr) - (n7) 

The basis {u n } n= _ oa with 

v n (x) = („!i„ (x) (118) 

is now given by 

Un (x) = Jn(27T|x|) _ ptf(x) „ _ (n9) 

The coefficients are related by 

a n = £ n a n ■ (120) 
Thus, for this example, the transformation matrix T is diagonal. The matrix elements R m n = 
E {<x m a* } and R mn = E {a m a* } are related by 
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3.4 Solving the Eigenvalue Problem 



The goal of this subsection is to approximate the eigenvalue problem for R by the eigenvalue 
problem of a finite-dimensional matrix. We shall replace the operators R and T in Equation (|103p 
by the corresponding truncated matrices and then use Equation f) 104(1 for these finite-dimensional 
matrices to find the eigenvalues. We shall state rigorous upper bounds for the approximation error 
of the corresponding eigenvalues. 

We follow the path indicated in [T5|, where the disk aperture of Example 2 has been treated. As 
shown above, the transformation matrix T is diagonal for that aperture. For the general case where 
T is not diagonal, the truncation problem is more involved. Moreover, it remains to be clarified 
how the truncation to finite-rank operators influences the diversity spectrum. In the following 
treatment, we shall answer these question. 



3.4.1 Truncating Errors of the Series 

As shown in [T3], the infinite series expansion of a (x) inside a disk is essentially finite-dimensional. 
The authors prove that the (absolute as well as the means square) truncation error that occurs by 
replacing the series (|115[) by the finite series 

N 

a(x;Af) = J2 a n v n (x) (122) 

n=-N 

decays exponentially with N. These results are based on estimates on Bessel functions that are 
studied in detail in Appendices I and II of |15) . The following theorem summarises the main results 
of those appendices. 



Theorem 1: Let r\ > and Njj = [enTil, where [.] denotes the ceiling function. Then, for 
N > Nd , the following two bounds hold uniformly in r for < r < r\ : 



J2 I J ™ (2^)1 < 0.2 exp [Nd - N] 



\n\>N 



Jl (27rr) < 0.01 exp [2 (N D - N)} 

\n\>N 



(123) 
(124) 



The first bound (|123[) controls the absolute truncation error when approximating a (x) by a (x; TV). 
Furthermore, it controls the error when p (x) given by the the series (|28[) is approximated by its 
truncated version 

N 

p N (x)= ~Sn^ n j n ^n(2nr) (125) 

n=-N 

This latter approximation is practically important because in general there is no closed analytical 
expression for p(x) available. Since the coefficients s n are bounded by (|27|) . the truncation error 
for r < ri and N > No = |"e7rri] can be estimated as follows: 



\p{x)-p N (x)\ = 



J2 S n e^ n j n J n (2nr) 

i\>N 
\n\>N 

0.2e ND ~ N 



(126) 



16 



The second bound (|124|) controls the means square (MS) truncation error of the series expansions 
for a(x). For Examples 1 and 2, the MS truncation error of the series (|108[) and (|115p can be 
estimated as follows: 















2' 


1 


E ® nvn 

\n\>N 




H 




\n\>N 


> 



= E {l«»l 2 }l^l 2 



\n\>N 

= E i&i a ( 127 ) 

\n\>N 

Here we have used v n = (, n u n , the orthonormality of the base u ni and the property 

E {|a n | 2 } = S = 1 

that follows from Equation ([50)) . In both examples, the coefficients £ n can be expressed by Bessel 
functions, see Equation (jllOp for the circle and Equation (|117[) for the disk. For the circle, |£„| = 
(27rr) holds, and we obtain the truncation error 







2 






E ® nV n 




> < E J »( 27rr ) 




\n\>N 




\n\>N 



< 0.01 exp [2 (N D - N)] 
with Njj = \cnr~\ . For the disk, we insert the expression (|117[) for £„ and obtain 







2^ 






E 




> = E iM a 




|n|>JV 


- 


\n\>N 



= I 1 ^rdr J « ( 2 ^) 

Jo \n\>N 

< 0.01 exp [2 (A D - TV)] / 7-j-rdr 

with A^d = [e7rri] . The integral equals one and we obtain the same bound as for the circle of the 
same radius: 







"1 




E ® nVn 


\ < 0.01 exp [2 (N D - N)] 




\n\>N 





This is just the statement of Theorem 2 in |15j . 

3.4.2 Approximations for Operators and Eigenvalues 

For the general case where the transformation matrix T is not diagonal, the formalism is more 
involved than in the two examples discussed above. Before going into details, we explain what 
problem has to be solved. 

Let R^r be the truncated version of the inifinite autocorrelation matrix R in which all elements 
Rmn with \m\ > N and |n| > N are set to zero. Heuristically, it is evident that R^r — > R in some 
sense as N — > oo. The appropriate type of convergence is specified and proven in Appendices A 
and B. From R^r — > R we then obtain convergence of the corresponding eigenvalue spectra. 
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We write Xi (Rat) and (R) for the eigenvalue number i of R^ and R, respectively. Assume that 
both eigenvalue spectra are numbered in descending order: 

Ai > A 2 > A 3 > ... (129) 

Recall that, since both operators are compact, they have discrete spectra. We shall have to prove 
that Xi (Rat) — > Xi (R) holds for N — > oo. Based on the estimates from [TS] as stated in Theorem 
1, we are able to prove the following 



Theorem 2: Let A C V where T> is a disk of radius r-y. Assume that S (a) is a piecewise 
continuous, bounded function over the interval [— tt,tt]. Then the inequality 

\Xi (R) - Aj (Rjv)| < 0.2 Pmax exp (N D - N) (130) 

holds for N > Nu and Njj = \eirri] and with a constant p max = max a£ [_ 7r j (27r<S (a)). 
The proof is given in Appendix B. 

Now let Tjv and Gat = T^Ta? be the truncated versions of T and G, respectively. Then Rat = 
TyvRTjy, and the approximate eigenvalues can be numerically calculated from 

cig (Rat) = cig (RGat) . (131) 

Examples of such numerically calculated eigenvalue spectra will be shown in the next section. 

According to Theorem 2, the eigenvalues of R can be approximated exponentially tight by the 
eigenvalues of Rat. One might expect that this has the consequence that the diversity measure 
<jj (R) given by Equation (|89p can be approximated tightly by 

^(Rjv) = 1 2 ■ (132) 
ll Rjv ll 2 

In fact, as shown in Appendix 3, the following tight approximation holds: 



Theorem 3: Let A C V where T> is a disk of radius r\. Assume that S (a) is a piecewise 
continuous, bounded function over the interval [— 7r, 7r]. Then the inequality 

<0Ap 2 max exp (N D -N) (133) 

holds for N > No and Nd = \enri \ and with a constant p max = max Q , e r_ 7r)7r i (2nS (a)). 

The statement of the theorem means that the inverse diversity measure (R) = ||R||2 can be 

approximated exponentially tight by the corresponding truncated quantity w _1 (Rat) = ||Rat|| 2 • 
These quantities are of primary relevance, because for the calculation of the slope Eo of a MIMO 
system according to Equation ((5|), these inverse diversity measures for the transmit and receive 
apertures have to be added. The approximation error of a numerically calculated diversity measure 
ui (Rat) itself can be estimated from the geometric series argument 

uj (R) = uj (Rat) — !— = ui (R w ) (l + e + O U 2 )) (134) 
1 — e 

with an error e = cu (Rat) ( ||Rjv|| 2 — W^-Wl ) ^ n at can be controlled by Equation (|133|> . 
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A 




Ay 



yx : Antenna 1 

y-2 : Antenna 2 

2/3 : Antenna 3 

z/4 : Antenna 4 



Az = vAt 

Figure 1: Four antennas mounted on a vehicle moving with velocity v. 



3.5 A Generalisation: M parallel Lines 

Let us consider M antennas mounted transverse across the roof of a vehicle that is moving into x- 
direction, see Figure [TJ Then the aperture A consists of M parallel horizontal lines. Time diversity 
provided by time interleaving corresponds to the x-direction, antenna diversity to the y-direction. 

Because this aperture is not a continuous curve, the results above can not applied directly. However, 
there is a simple argument how they can be adopted to this important case. 

The problem can be solved by constructing a Hilbert space T-L on that aperture that consists of M 
straight lines. The construction is just a slightly modified direct sum of the M Hilbert spaces 



1, M corresponding to these M straigth lines. Let (ft m , ift m £ TL m . Their scalar product 



hif, 



is denoted as (cft,ip) m . The vectors (ft, ift £ % are then defined as M-tuples (ft = (</>i, 
tp = ...,ipM) T that, for convenience, may be written as columns. Their scalar product is 
defined by 



M 

-Y 

m—l 



{(ftmAra) 



(135) 



Let r be the common priviledged parameter of all the lines. Then the scalar product can be written 

as 



M 

-Y 

m—l 



C(r)V(r)dr. 



(136) 



The fading process itself is a column vector a = (ai, aj/) ■ The corresponding autocorrelation 
operator R will then be given according to the general definition (|75|) as 



,R0) = E {{<f>, a) (a,V)} 



(137) 



The corresponding integral kernel now becomes a tensor JZ mm i (t, r'), m,m' £ {l,...,Af}. This 
kernel is a continuous matrix- valued function in t and r'. As a consequence of that continuity, all 
the results above remain valid and our method can be applied to this case. 



4 Numerically calculated diversity spectra 
4.1 How to Calculate Diversity Spectra 

The task now is to calculate the eigenvalue spectrum according to Equation (|13ip . It depends on 
the geometry of A and on the Fourier coefficients s n of the PAS S (a). The elements of Gjv are 
given by G mn = (v m ,v n ). This scalar product is an integral over A and has to be calculated by 
numerical quadrature methods for all m, n with |m| , \n\ < N. The matrix elements of R are given 
by R mn — 27rs m _„. The truncation number N governs the accuracy of the eigenvalues according to 
Theorem 2. For the choise N = Nn + 10, e.g., the exponential is approximately 5 • 1CP 5 , which can 
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be regarded as sufficiently accurate. It is noteworthy that the constant 0.2 • max a6 [_ Ijr ] (2irS (a)) 
in that theorem is typically in the order of onfl The number No = [Wri] grows with the size of 
the aperture A. If, e.g., it is inside a disk with radius of one wavelength, No = \ctt~] = 9. Thus, 
for typical aperture restricted to the size of a few wavelengths, the matrices do not exceed the size 
of a few times ten. 



4.2 The Power Azimuth Spectra 

We consider two PAS prototypes that are quite popular in the literature: The uniform PAS and 
the von-Mises PAS. 



Uniform PAS: By this we mean a PAS that is constant over an opening angle A. The normalised 
and centered uniform PAS of opening angle A is given by 

S (a; A) = irect (|) . (138) 
The corresponding Fourier coefficients can be calculated as 

1 -sinc ( n— ) . (139) 



2?r V 2tt 



von-Mises PAS: The von-Mises distribution is a circular distribution with a shape quite similar 
to a wrapped normal distribution |24| . Its shape is controlled by a parameter k > 0, whose inverse 
is approximately equal to the variance of the corresponding normal distribution. The centered 
PAS is given by 

S(a; K ) = —^—e—\ (140) 
27rl (k) 

where I n (k) denotes the modified Bessel function. The Fourier coefficients 

In (k) 



2ttI (k) 

can be obtain from Formula 9.6.34 in |25) . 



(141) 



The formulas above apply to a PAS that is centered around a = 0. The versions centered around 
a = o;o are obtained by the shift 

S(a) ^S(a-ao) . (142) 
The corresponding Fourier coefficients must then be modified according to 

Sri. ' ^ ^rj.e ^ 



4.3 Uniform Circular Arrays 

A uniform circular array (UCA) is a configuration where the antennas are mounted uniformly on 
a circle. In the limit of a dense array, this is just the situation of Example 1 discussed above. 
Figures [5] and [3] show the numerically calculated eigenvalue spectra for a ULA of radius r = 1 (i. 
e. one wavelenght) for the uniform and the von-Mises PAS, respectively. The parameters A and k 
are adjusted manually in such a way that the corresponding values for w are nearly identical. The 
eigenvalue spectra are very similar, but, however, they are not identical. Obviously, the multipath 
richness decreases significantly for small opening angles of the PAS, i.e. for small values of A and 
for high values of n. 

Figures 2] and [5] show ui for small radii up to r = 0.5 for the same parameters of A and k. The 
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-O- A=45 
-B- A= 90' 




Eigenvalue index 



Figure 2: Diversity spectra for a dense UCA of radius r = 1 for a uniform PAS with different 
opening angles A. 




Eigenvalue index 



Figure 3: Diversity spectra for a dense UCA of radius r = 1 for a von-Mises PAS with different 
parameters k. 
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Figure 4: The diversity measure w as a function of the radius for a dense UCA and a uniform PAS 
with different opening angles A. 




Figure 5: The diversity measure w as a function of the radius for a dense UCA and a von-Mises 
PAS with different parameters k. 
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Figure 6: The diversity measure w as a function of the size Ax for a dense ULA and a uniform 
PAS with different opening angles A. The main direction of arrival is ao = 0, i.e. parallel to the 
array. 

curves are nearly identical in both figures. The green curve of the isotropic PAS (A = 360° or 
k = 0) shows a noticable oscillatory behaviour, that can also be observed for other values that are 
close to the isotropic case. 

4.4 Uniform Linear Arrays 

For a straigth line which corresponds to a uniform linear array (ULA). such an oscillatory behaviour 
could not be observed. Figure [5] shows ui for ULA sizes up to one wavelenght and the uniform PAS. 
The direction of the ULA is parallel to the mean angle of arrial (ao = 0). This is the worst case, 
and for small values of A (blue and red curve), there is practically no diversity gain. Increasing 
the size of the ULA up to several wavelengths does not really help because the steepness of the 
curves is very poor. There is a strong dependency on the mean angle of arrival ao- Figure [7| shows 
w as a function of «o for Ax = 10. It can be seen from this figure that even for this large array, the 
diversity measure falls down close to one for small opening angles (blue and red curve). This means 
that for a narrowly focussed signal, even a huge antenna array parallel to the direction of arrival 
can hardly provide any diversity gain. Such narrowly focussed signal may be not very likely at the 
mobile station, but they will frequently occur at the base station site. Therefore it is important at 
that site not to fix all antennas along one line, but to distribute them in two dimensions by using, 
e.g., a circular array or any other two-dimensional configuration. Figure [5] shows u> as a function of 
ao for a configuration of four parallel ULAs of length Ax = 1 and a separation of the outer ULAs 
given by Ay = 1. In that scenario, the diversity measure is very robust against rotations. As 
discussed above (see Figure [TJ), by space-time scaling Ax = vAt, this situation can also be related 

6 For the isotropical PAS S (a) = (2n)~ 1 , as an example, it has the value 0.2. For a uniform distribution 
S (a) = A - red (a/ A) over an opening angle A, it has the value 0.2 • 2n/A, which does not exceed the order of 
one unless the opening angle is extremely narrow. 
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Figure 7: The diversity measure w as a function of the main direction of arrival a® for a dense 
ULA and a uniform PAS with different opening angles A. The size of the array is Ax = 10. 



4<>- 



A= 45 
A=90° 
A=180° 
A=360° 




Figure 8: The diversity measure w as a function of the main direction of arrival ao for four parallel 
and equidistant lines and a uniform PAS with different opening angles A. The length of the lines 
is Ax = 1 and the separation of the outer lines is Ay = 1 . 
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Figure 9: The diversity measure us as a function of the number of antennas for a UCA of radius 
r = 1 and a uniform PAS with different opening angles A. 

to the scenario of a vehicle moving in x-direction with velocity v and a 4-antenna communication 
system with time interleaving over a interval At. 

4.5 Discrete Arrays 

The focus of this paper is to analyse the diversity spectra of continuous (dense) arrays, i.e. of 
contiuous regions in the plane. It is interesting to compare this with the discrete case where a 
finite number L of antennas is mounted on a finite geometrical base, i.e. a circle for the UCA or 
a line for the ULA. In [13], a detailed analysis has been presented how the diversity measure io 
depends on the array geometry. That analysis, however, is restricted to ideal case of an isotropic 
PAS where the autocorrelation function is given by po (x) ~ Jo(27rr). In the general case, the 
ACF is given by the expansion (|28|) . According to Theorem 1, this can be approximated by the 
truncated series 

N 

PN (x) = ]T s n ^ n j n J n (2nr) . (143) 

n=-N 

The appropriate order N can be chosen as N = No + 10 according to the bounds obtained in 
the last section. The diversity measure w according to Equation has be calculated from the 
matrix R with elements 

R lk = p (xj - x fe ) w p N (xj - x fe ) , (144) 

where x^, i = 1,...,L are the antenna positions. We note that Equation (|86p does not require 
to solve the eigenvalue problem for R. Only the traces of the matrices R and RtR need to be 
calculated. 

Figure shows the depency of w on the number L of antennas for a UCA of radius r = 1 and a 
uniform PAS of different opening angles A. The margenta dash-dotted line corresponds to oj = L 
(i.e. the diversity measure of L uncorrelated antennas). The horizontal dashed lines mark the u>- 
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Figure 10: The diversity measure uj as a function of the number of antennas for a ULA of size 
Ax = 2 and a uniform PAS with different opening angles A. The signal comes from broadside 
(oo = tt/2). 



values for the continuous UCA and may be interpreted as the asymptotic limit L — > oo. It can be 
seen from the figure that for large L, the diversity measure runs into a saturation given by that 
continuous case asymptote. It is interesting to note that uj does not increase monotonically with 
L. This interesting fact has been already observed in [T3] for the isotropic case (i.e. A = 36O°0. 
In the same paper, it has also been pointed out that for a ULA, the maximal value of uj is typically 
reached for a relatively small number of antennas, and the limiting value for L — > oo lies below 
that value. Figure [TO] shows the situation for a ULA of size Ax — 2, ao = 90° and the same PAS 
parameters as above. For A = 360°, the maximal value of uj is reached at the number of L = 6 
antenna. For A = 90° and A = 45°, this number is even lower. The black curve for A = 180° 
hidden by the green curve corresponding to A = 360° because of the mirror symmetry of that 
geometrical configuration as already discussed in the context of the Doppler effect. Figure [TO] (or 
similar figures for other values of Ax ) provides practical design guidance for the design of linear 
arrays. The number of antennas should be chosen in such a way that uj is not only maximal for a 
special case like A = 360° but should stay robust against a change of the PAS. 



5 Discussion and Conclusions 

The diversity corresponding to an aperture A is characterised by the diversity spectrum {A}. , of 
eigenvalues of the spatial autocorrelation operator R that acts as an integral kernel on the Hilbert 
space of square integrable functions defined on that aperture. The diversity spectrum depends on 
the geometry of the aperture and the power azimuth spectrum that characterises the statistics of 
the spatial fading. We have shown a method for the numerical calculation solving the eigenvalue 
problem for an finite-dimension approximation of R. The method provides rigorous bounds of the 
approximation error. A simpler, but less comprehensive characterisation compared to the whole 

7 The green curve of Figure [9] is identical to the black curve of Figure 4 in that paper. 
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diversity spectrum eig (R) = {A} i=1 is a single number named the diversity measure w (R). It can 
(roughly) be interpreted as the effective number of equivent independent diversity branches. Its 
calculation is in the same framework, but it is simpler because no eigenvalue problem has to be 
solved. The dependency of to (R) on the aperture geometry, on the main angle of incidence, and 
on the anisotropy of the PAS gives useful hints for the design of antenna arrays. For a vehicle 
of constant velocity moving through the spatial fading pattern, the time dependent fading can be 
obtained by a simple scaling, and a time interval corresponds to a spatial line segment. Therefore, 
the time interleaving depth for a coded system can be assigned to a line segment, and the diversity 
of that segment corresponds to the code diversity that can be exploited within the time interleaving 
depth. 



A Some Definitions and Facts about linear Operators 

Let A be a linear operator on H. Its (hermitian) adjoint, , is defined by 

(0,Aty) = (A^). (145) 

If A = A^ holds, the operator is called self-adjoint. In this paper, we need only to consider 
bounded operators. For these operator, self-adjoint is the same as the (in general weaker) property 
symmetric or hermitian. 

A limit R^ — > R for operators may be defined in different meanings. In this paper, we need three 
types of norms for an operator A on a Hilbert space H. We briefly recall their definitions and their 
most important properties and refer to j20j for further details. 

1. The {uniform) operator norm ||A|| is defined by 

\\H = sup ■ (146) 

11011=1 

Operators with ||A|| < oo are called bounded operators. 

2. The Hilbert- Schmidt norm ||A|| 2 is defined as 



||A|| 2 ^tr(AtA), (147) 

where tr(-) denotes the trace of an operator. For the special case of finite-dimensional 
matrices, the Hilbert-Schmidt norm is called Frobenius norm. One can show that 



|aii 2 = ./Ei^i 2 ( 148 ) 



/ k 



is an equivalent definition. In that equation, Aik = {ui,Uk) are the matrix elements in 
an arbitrary orthonormal basis Ui. Operators with ||A|| 2 < oo are called Hilbert-Schmidt 
operators. The Hilbert-Schmidt norm can also be written as 



|A|| 2 = 



,E A ?> ( 149 ) 



where the Af are the eigenvalues of A^A. In case that A is self-adjoint, the real numbers Aj 
are just the eigenvalues of A. 

3. The trace norm HA^ for a self- adjoint operator is given by 

oo 

\\A\\ 1= J2\^\, (150) 

where Ai are the eigenvalues of A and a purely discrete spectrum of A has been assumed. 
Operators with HAj^ < oo are called trace class operators. 
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In our paper, we make use of some properties of compact operators. For the exact mathematical 
definition of compactness, we refer to [20] ■ For our purposes, it is sufficient to recall that every 
Hilbert-Schmidt operator is compact (see Theorem VI. 22 in |20|). 

From the above definitions between norms, the relation 

HAll^llAtAll^llAH? (151) 

follows immediately. Furthermore, there is the following hierarchical relation holds between the 
norms defined above ([5D], Theorem VI. 22): 

llAHIIAll^HAll! (152) 

The following inequalities also hold 

IICABU, < IICUIIAIIJBII , (153) 
||CAB|| 2 <||C||||A|| 2 ||B|| , (154) 

provided that all the norms exist. 



B Proof of Theorem 2 

In this appendix, we proof the bound (|127[) for the approximation error of the eigenvalues as stated 
in Theorem 2. This bound is based on bounds for approximation errors for operators. 

First of all, we shall proof that the operator R with matrix elements R mn given by Equation (1301) 
is bounded. 



Lemma 1: Under the assumption that S (a) is a piecewise continuous, bounded function over 
the interval [— 7r, 7r], the operator R with matrix elements R mn given by Equation pop is bounded 
with operator norm given by 

llRll = 2tt max S (a) = p max 

a£ [— 7r,7r] 



Proof: Since we have assumed that S (a) is a piecewise continuous, bounded function over the 
interval [— 7r,7r], it defines a bounded multiplication operator 

S : X (a) ^ Y (a) = S (a) X (a) 

on C 2 ([— 7r, 7r] , |^). Its operator norm is given by 

||S||= max 5(a). 

The Hilbert space C 2 ([— 7r, 7r] , 4 s ) is isomorphic to the Hilbert space £ 2 (Z) of square summable 
sequences with indices running over Z. The vectors X (a) , Y (a) £ C 2 ([— 7r, 7r] , 4 s ) correspond to 
the sequences {x n } neZ , {y n } n£ z ^ ^ 2 (^) °f their Fourier coefficients. The action of the operator 
S in this discrete representation is given by the convolution 

oo 
m— — oo 

where s„ are the Fourier coefficients of S (a). Comparing with Equation (|3"0|) we see that R = 27rS 
holds in this represention. Since the operator norm does not depend on the representation, we 
conclude that ||R|| = 2tt ||S||, which completes the proof. ■ 

For two compact and self-adjoint operators, the difference of the corresponding eigenvalues is 
bounded by the operator norm of their difference. We cite the following fact from the textbook of 
Riesz and Nagy [22], § 95: 
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Lemma 2: Let A and B be self-adjoint compact operators and define 

C = A + B. 

Assume that the eigenvalues Xi (A) of A and Xi (C) of C are ordered according to (|129p . Then 
their difference is bounded by 

IMC) -A, (A)| < ||B|| 

for all i. 



Corolarry to Lemma 2: Setting C = R, A = R^v, and B = R — R^r leads to the inequality 

|A i (R)-A i (R JV )| < ||R-Riv]| (155) 
The goal now is to show that ||R — Rat|| — > holds for N — > oo. 

In the following, we assume that A C V, where I? is a disk of radius n. We define Nrj = \cnri~\ 
and assume that N > Nrj- 

The transformation matrix T is defined by its elements T mn = (u m ,v n ). The Gram matrix 
G = T^T has matrix elements G mn = (v m , v n ). The squared Hilbert-Schmidt norm of T is given 
by 

oo 

||T|| 2 = 11^^ = 11011!= £ ^v n ). (156) 

n— — oo 

Using Equation (|97p we may write 

(v n ,V„) = / \v„ (x)| 2 d/i(x) . 
J A 

J I 2 l (27T|x|)d Ai (x) 

A 

The integral is over the array A with dfi (x) given by Equation ([Ml) or (|63|) for the respective cases. 
Inserting the above expression into Equation (j!56[) leads to 



ITII 2 



The properties 



and 



then yield 



/ £ J 2 (2 7 r|x|)d M (x) 

n— — oo 
oo 
n— — oo 

/ d/i(x) = l 



|T|| 2 = 1. (157) 



We now define a truncated transformation matrix Tjy that has the same matrix elements as T 
inside the middle square defined by |m| < N and \n\ < N. All other matrix elements are set to 
zero. The remainder matrix T — Tjv has only non- vanishing matrix elements outside that middle 
square. Its squared Hilbert-Schmidt norm is given by 

\\T-T N f 2 = tr [(T - T N ) f (T - T N ) 

= tr (T+T + Tj/IV - T^T N - Tj/T 



tr (T*T) + tr [ T\ t T n ) - tr (T^T N ) - tr ( TlT 



tr (T^T) -tr (Tt r T 



tr ( T*T - Tt r T 



AT J- AT 
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Here we have used the fact that T^TV = T^TV and thus tr (T^Tjv) = trfTjyTVj. The 

truncated Gram matrix Gat = T^Tat has the same matrix elements as G inside the inner square 
and zero elements outside. Thus, the expression derived above can be written as 



HT-Tjvlla = tr(G-Gjv) 

\n\>N 



f ]T 4(27T|x|)d M (x) 

JA \n\>N 



From Theorem 1 we get the bound 



Taking the square root yields 



From 



|T - T^v H2 < 0.01 exp [2 (N D - TV)] 



IT-TjvIIj, < 0.1 exp (N D - N) 



(158) 



\n\<N 

and Equations (|156[) and (| 157[) we conclude that 

\\T N \\ 2 < ||T|| 2 = 1 
holds. We define the truncated correlation matrix 

t 



Rjv — TatRT^y 



(159) 



(160) 



The Hilbert-Schmidt norm approximation error relative to R = TRT^ can now be estimated as 
follows: 



IR-R 



N\\ 2 



TRT+ - TArRTjv 



TR (T - Tat) 1 + (T - Tat) RT f 



N 



< 



TR (T T 



(T-TAr)RTj; 



N 



< ||T|| 2 ||R|| ||T - T W !| 2 + ||T - TaH| 2 ||R|| \\Tn\\ 2 

< 2||R|| ||T — Tjv|| 2 

From the bound (j!58[) we conclude 

||R - Rjv|| 2 < 0.2 ■ ||R|| exp (N D - N) . 
The norm hierarchy (|152p yields 

||R — R-jvII < 0.2 • ||R|| exp (N D - N) . 
From the Corollary to Lemma 2 we obtain 

\Xi (R) - A, (11^)1 < 0.2 ■ R exp (N D - N) . 

Lemma 1 then leads to Theorem 2. 



(161) 



(162) 
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C Proof of Theorem 3 



In this appendix, we prove the bound (| 133[) for the approximation error of the squared Hilbert 
Schmidt norms as stated in Theorem 3. From R = TRT 1 and R N = T^RT^ and Inequality 
(11591) we conclude 

||R|| 2 < ||R|| and ||Rjv|| 2 < ||R-|| ■ 



Then 



|R|| 2 - ilRwils 



The bound (flFTj) yields 



IRIU- IIR 



Jv || 2 



|R|| 2 + ||Rjv|| 2 ) 



< 2 R 



|R|| 2 - ||Rjvil 2 



|R|| 2 -||Riv|| 2 



< 0.4||R|| 2 exp(^ I) - N) 



With Lemma 1 we then obtain Theorem 3. 
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